Magnetism and homogenization of micro-resonators 



Arrays of cylindrical metal micro-resonators embedded in a dielectric matrix were pro- 
posed by Pendry, et. al. , [T7J as a means of creating a microscopic structure that exhibits 
strong bulk magnetic behavior at frequencies not realized in nature. This behavior 
arises for f/-polarized fields in the quasi-static regime, in which the scale of the micro- 
structure is much smaller than the free-space wavelength of the fields. We carry out 
both formal and rigorous two-scale homogenization analyses, paying special attention 
to the appropriate method of averaging, which does not involve the usual cell averages. 
We show that the effective magnetic and dielectric coefficients obtained by means of 
such averaging characterize a bulk medium that, to leading order, produces the same 
scattering data as the micro-structured composite. 

Keywords: homogenization; meta-material; micro-resonator; magnetism; quasi-static. 

1 Introduction 

Within the field of artificial materials, there is presently intense activity in the area of 
creating "metamaterials" with negative bulk dielectric or magnetic response. Materials with 
dielectric or magnetic coefficients that are either simultaneously negative or of opposite sign 
offer a rich variety of interesting and useful phenomena. As nature provides us with materials 
that exhibit negative response only at rather restrictive frequencies, one of the aims of this 
field is to extend the selection of frequencies by creating microscopic structures that have 
resonant response when natural materials tend to be unresponsive. The field received a 
jump start with the introduction of model structures of thin metallic wires for creating 
electric resonance [16] and ring-type structures for creating magnetic resonance [17], which 
were proposed by Pendry, Holden, Robbins, and Stewart in the late 1990s. Combinations 
of these effects were investigated by Smith, et. al. [21] and many others, to create "left- 
handed" materials, possessing a negative index of refraction. More recently, Pendry [T5] 
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Abstract 



proposed the creation of negative refraction by composites in which one of the components 
is chiral, and the homogenization of such structures has been investigated in [llj. There 
seems to be considerable debate and some confusion concerning the definition and meaning 
of bulk effective electromagnetic coefficients in this setting; moreover, rigorous mathematical 
treatment of the subject is still in its early stages. 

Our intention with this work is to help clarify the meaning of the effective bulk dielectric 
permittivity and magnetic permeability in the quasi-static limit, in which the scale of the 
micro-structure is small compared to the free-space wavelength of the fields. We work only 
with a two-dimensional model of ring-type resonators, for which magnetism is the dominant 
effect. The mathematical context of our study is periodic homogenization. The relations 
between the D and E fields and between the B and H fields of the individual components 
of the micro-structured composite material, 

D = sE, B = /iH, (1) 

give rise to bulk coefficients e* and /x* that govern certain average fields on the macroscopic 
level: 

D av = e*E av , £? av = /x*if av . (2) 

What is noteworthy in the homogenization of micro-resonators is that these average fields 
are not to be understood in the standard way as micro-cell averages. Indeed, as Pendry, 
et. ai, p2] observed, even if /x = /xo for all components, we still obtain a nontrivial magnetic 
response, /x* 7^ /Jo- 

The crucial ingredient for emergence of magnetic behavior is the presence of a component 
with extreme physical properties. In fact, the rings in our resonators must possess high 
conductivity or internal capacitance tending to infinity as the inverse of characteristic length 
of the micro- structure. 

For broad discussions of electromagnetic materials with negative coefficients, one may 
consult [18], [23], or [19], for example. 



1.1 Magnetism from micro-resonators 

In our model, the micro-resonators are represented by infinitely long rods with conducting 
surfaces (Figure [T]). The fields are harmonic (with frequency uS) and magnetically polarized. 
The magnetic field, denoted by the scalar h(x\, x 2 ), is directed parallel to the rods, while the 
electric field E(x\,X2) lies in the plane perpendicular to the rods. The electric field induces 
a current j on the surfaces of the resonators; it is related to the tangential component 



E(x\,X2) through a complex (more on this in section 1.3) surface conductivity cr s : j = a s E-t. 
The current, in turn, effects a discontinuity in the magnetic field, h x — h e = j, where h e 
and hi denote the values of h exterior and interior to the micro-resonators. The Maxwell 
system of partial differential equations reduces to the following system for h and E (where 
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kx V = (—d/dx 2 ,d/dxi)): 

off the surfaces of the micro-resonators, (3) 



V • E — iu/ih = 0, 
V x h - iuoeE = 0, 

h e — hi + a s E-t = 0, with E-t continuous on the surfaces, (4) 



with the convention that the unit tangent vector is directed in the counter-clockwise sense. 

The quasi-static limit amounts to fixing the frequency and allowing the period of the 
micro-structure to tend to zero. This is to be contrasted with work of Sievenpiper, et. ai, 
[20] for example, who devise capacitative structures for the manipulation of photonic spectral 
gaps, a phenomenon that is pronounced when the wavelength and period are comparable. 
Starting from ([3]-[4]), we shall derive a system of Maxwell equations governing suitably defined 
macroscopic fields. Because the currents around the resonators flow in microscopic loops, 
they do not appear as currents in the homogenized equations. Instead, their bulk effect is 
manifest through the effective magnetic coefficient //*, which is complex (even if e and /x are 
real) because it incorporates the effect of loss due to the currents. Our homogenized system 
is 

V^av-iw/i*/^ = 0, 

V ± h° c -iue*E ! , v = 0. U 

It is important that the equations for the bulk fields are posed in terms of the exterior value 
of the magnetic field, which is denoted by h®, not in terms of its cell average. 
At the risk of redundancy, we emphasize the following two key points: 

1. The characteristic of the micro-structure that is crucial for the emergence of magnetic 
response from nonmagnetic components is that one of the material properties, the 
surface conductivity, is extreme. More precisely, it scales inversely with the microscopic 
length scale. 

2. In the homogenized Maxwell system, the macroscopic H field is not the cell average 
of the field but rather the value exterior to the resonator. The macroscopic B field, 
however, is the usual micro-cell average. 

In fact, not only are the H and B fields averaged over different parts of a unit cell, but 
so are the E and D fields. This means that, just as we have discussed for the magnetic 
coefficient, even if e — E\ in all components of the unit cell, the effective dielectric coefficient 
e* will typically be different from e\. This feature distinguishes the homogenization of micro- 
resonators from the more "standard" homogenization of composites with perfectly bonded 
interfaces and material properties that are not extreme, in which all macroscopic fields are 
understood as micro-cell averages. 

The second feature is already present in the problem of homogenization porous media, 
in which the cell average is taken outside the holes [8] . More recently, both of these features 
appeared in the work of Bouchitte and Felbacq (U El QUI E] , who demonstrated the emergence 
of bulk magnetic behavior from nonmagnetic materials in a somewhat different but related 
problem. The extreme property in their setting is the dielectric coefficient inside a periodic 
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inclusion, which tends to infinity as the inverse area of a micro-cell. The magnetic field in 
the matrix, exterior to the inclusions, has vanishing fine-scale variation and appears as the 
macroscopic field in the homogenized equations. This exterior value drives the fine-scale 
oscillations in the interior, whose resonant frequencies produce extreme magnetic behavior. 
The point in their problem as well as ours is that the effective equation involves only the H 
field exterior to the inclusion (or resonator), but the B field must be averaged over the entire 
micro-cell. A similar problem in which fiber arrays have extreme conducting properties in 
the plane perpendicular the fibers but not in the direction of the fibers has been investigated 
rigorously by Cherednichenko, Smyshlyaev, and Zhikov [7j. Here again, the field in the 
matrix appears in the effective equations, which possess the additional feature of spatial 
non-locality in the direction of the fibers. These problems are to be contrasted with the case 
of a small-volume-fraction array of conducting metallic fibers of finite length, treated by 
Bouchitte and Felbacq |Bj. In that setting, the conductivity is extreme but the field averages 
are nevertheless taken in the usual way. 

We shall discuss the homogenized system ^ — and our scheme for defining "averaged 
fields" and "effective properties" — further in Section|2j and we offer a systematic justification 



in Section 4.4 following the formal asymptotic analysis. The main points are these: 



1. The effective coefficients e* and fi* describe the limiting behavior of the scattering 
problem. 

2. Taking the value of H exterior to the resonators as the macroscopic field is the only 
choice that preserves the Maxwell-type structure of the effective system of PDEs (see 



the comments following equation (19)). 



3. The averaging scheme is consistent with the treatment of the E, D, H, and B fields 



as differential forms (section 4.4). 



Point (1) means that the field scattered by a micro-structured object when illuminated 
by a plane wave should tend to the field scattered by an object of the same shape consisting 
of a material possessing the bulk coefficients. In particular, the reflection and transmission 
coefficients associated to scattering by a micro-structured slab should approach those for the 
homogenized slab. This is the model we have chosen for our analysis. It is important to 
keep in mind that, in any scattering problem, we assume that no resonator is cut or exposed 
to the air, so that the H field in the matrix dielectric surrounding the resonators connects 
continuously with the H field in the air. In other words, the air-composite interface cuts 
only through the matrix material. Although we do not treat boundary-value problems, it is 
evident by the same reasoning that the same effective equations remain valid for problems 
in which the boundary is exposed only to the matrix. 



1.2 Scaling of fields in the quasi-static limit 

Let us take a more careful look at our particular scalings. 

As one expects, the variation of h at the scale of the micro-structure, which we denote 
by f], vanishes in the matrix as this fine scale tends to zero (the quasi-static limit); the 
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microscopic variation of h is due solely to its discontinuities at the current-carrying surfaces 
of the resonators. (On the other hand, the electric field E will have non-vanishing micro- 
periodic oscillations in the matrix, even if e and fi are constant.) In the quasi-static limit, 
therefore, the jump in h, which we call the current j, is constant around a single micro- 
resonator. 

We have emphasized the point that interesting magnetic behavior arises only when the 
micro-resonators are highly conducting. This means that the surface conductivity cr s tends 
to infinity as the size rj of a period cell tends to zero. The reason is described by Pendry, 
et.al, [T7J and is borne out by our homogenization analysis. The electromotive force (EMF) 
around a single resonator is given by the line integral of the electric field around its surface. 
If the interior of one resonator occupies a region G v , then, using the first of the Maxwell 
equations pi), we can write the EMF in two ways: 



As h is of order 1 in 77, E • t is forced to be of order 7] on the surface of the resonators. 
Now, using the relation h e — h x + a s E ■ t = 0, we observe that, in order that nontrivial 
behavior emerge at the microscopic level, we should use a highly conducting material, namely, 
a s ~ v(~ x . In summary, the relevant scalings in our analysis are the following. 



1.3 The model for micro- resonators 

If the resonator consists of a solid metal cylinder, the connection between the current and 
the electromagnetic fields is accomplished by a simple constitutive law relating j to E ■ t 
through a real scalar a s , the surface conductivity. In this case, the resonator acts purely as 
an inductor. If the solid cylinder is replaced by a uniform solid metal ring, effectively the 
same situation persists, as there is no uneven distribution of charge on the two surfaces of 
the ring that would give rise to capacitative effects. 

More elaborate resonators allow for capacitative effects by forcing nonuniform build-up 
of charge around the ring. These are known as split-ring resonators (SRR) (see [17] or |19j . 
for example), composite ring structures consisting in part of dielectric material and in part 
of one or more incomplete metal rings. As the temporal distribution of charge is out of phase 
with that of the current, so also are the inductive and capacitative contributions to the jump 
in h. 

It is not our intention to provide a rigorous account of the inductive and capacitative 




(6) 



(7) 
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effects of SRRs in the quasi-static limit]]] Rather, in this work we are content to observe 
that in the literature, the electromotive force (EMF) around the ring arises from integrating 
two quantities around it, both related to the current j: 

1. An inductive component pj, where the real quantity p is the effective resistance of the 
metal in the composite ring. 

2. A capacitative component j/(—iuC), where the real quantity C is the capacitance 
arising from splits in the rings or the proximity of two closely stacked metal rings 
within the compound ring. 

Now, using the representation EMF = J dG E-tds, we model a general resonator by a 
single closed loop (a cylindrical surface in the three-dimensional realization) on which the 
inductive and capacitative effects are manifest through a single phenomenological complex 
constitutive law j = <J S E4, where a s = + ia s 2 . This amounts to defining a singular current 
a^E ■ t and a singular D field of strength —{a^/uj)E ■ t around the resonator. 

We now demonstrate the nature of the correspondence between our formula for the 
effective magnetic permeability and that obtained by Pendry, et. al. First, consider the case 
that p = po in all components and that the resonator is represented by a simple metal 
cylinder of (nondimensional) radius R < 1 in relation to a scaled unit cell (Fig. [2j right). 
The actual radius is r = rjLR, where L is an arbitrary fixed length corresponding to 1 in the 
macroscopic variable x. Keeping in mind that the conductivity should tend to infinity with 
rj" 1 , we set it equal to 

<h = , (8) 

VP 



where p is a fixed real constant. From equations (59), (61), and (68), we obtain 

2ip 



P = Po 



1 — ixR 1 + 



ujRpo 



(9) 



which is the formula (13) in [T7]jj 

Now, in order to incorporate capacitance into a SRR as in Figure [2j we take a phe- 
nomenological step and observe that, in order that our model produce formula (20) of [17J, 
the complex conductivity must now be set equal to 

1 \ ( 3A \ _1 
a s = -{p + iry 1 = - (p + i— , (10) 



1 The task of providing such an account, i.e., justifying equation (10) below starting from the Maxwell 
equations, remains in our view an open problem worthy of analysis. 

2 To make the connection to formula (13) fj, e g = 1 - ^ (l+ and (17) or (20) fi cS = 1 - 

( x _|_ — 3d — j m rj2] one mus t relate the notation in that paper to our notation in the 

following way: a r/L, r h- > r, d i— > d, /logoff t— > ^* i a l— * V^P- Note that the symbol a in |17j denotes the 
resistance of the metal, which is the reciprocal of conductivity. 
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where d = rjLA is the small distance between the outer and inner shells making up the 
circumference of the resonator. Here, A is the nondimensionalized distance in a rescaled 
unit cell (Fig. [2j left). The coefficient €q appears in 10 because it is assumed that e = €q in 
the dielectric between the inner and outer rings. The formula for /i* becomes 



1 - nR 2 1 + 



(ip- 



In order that the imaginary part of a s tend to infinity as rj 1 , we must have 

A~r/ 2 . (12) 

This means that the actual distance d = t]LA between shells would have to decrease as the 
cube of the cell length in order for the complex part of a s to have an effect in formula (|9j in 
the limit as rj — > 0. Of course, if the dielectric coefficient is allowed to tend to zero at some 
rate, the extreme rate of convergence of A to zero can be relaxed. 



We observe from (11) that if /xq is real and positive then /i* has positive imaginary part 



(provided p > 0). The assumption used above that pi = fj,Q in all components was merely 
a simplification; even if \i varies in space, a similar calculation of \i* is straightforward. It 
reveals that, if \i is real and positive throughout the structure, then 



Im(/i*) > 0. 



(13) 



We shall need this fact in Section HI 

We note that our homogenization-based analysis might not always be adequate for model- 
ing the behavior of a specific device. In particular, it is difficult to fabricate micro-resonators 
with shells that are extremely close to each other; moreover, the regime of interest is often 
not quasi-static, e.g., the wavelength may be only several times the length of a unit cell. 



1.4 Our approach 

The effective tensors e* and /i* and the unit-cell problem that determines them, together 
with the first-order correction to the leading order H field, are obtainable through formal 
asymptotic analysis, which we carry out because we feel it provides a clear intuitive point of 
view. We then prove the main results rigorously using the method of two-scale convergence 
(see Allaire [1] for a systematic treatment). 

The main ideas of the rigorous two-scale arguments are these. The weak form of the 
scattering problem, Problem [2] in Section [5]), is posed within the ^-dependent space 

H\ni)®H\^l (14) 

where Q% and are the domains exterior to and interior to the micro-resonators. The 
boundary term, involving the conductivity, brings the jump discontinuity of the H field into 
the equation. Notice that, if the conductivity remained bounded as rj — > 0, then this term 
would be of order ?7 _1 and the jump would disappear in the limit. Thus, as we have pointed 
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Figure 1: The strip S, consisting of one period in %2 of the micro-structured slab and 
surrounding air. The segments r± are artificial boundaries used in the weak formulation of 
the scattering problem. We define also f2 c = Vl \ f2 ; and fi 0c = \ ^i- 




Figure 2: Left: An example of a split-ring resonator. The current flows in the direction 
of the arrows. The splits in the rings cause the charge to be nonuniform, creating a large 
.D-field between the rings. This allows the current to follow a complete circuit by passing 
from one ring to the other through a capacitative gap. Right: The microscopic unit cell 
Q coordinatized by the variable y = x/rj. The idealized metal-dielectric micro-resonator 
consists of the boundary dG of a simply connected domain G. 
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out, we take the conductivity to be of order 77" 1 so that the boundary integrals remain on 
the same order as the the area integrals. 

We then show that the scattering problem always has a solution h v for each r]. As we do 
not have a priori bounds on the solutions, we first scale them and obtain the two-scale limits 
of the scaled solutions as well as those of their gradients. The uniqueness of the solution of 
the homogenized system governing these scaled solutions then allows us to obtain a postiori 
uniform bounds on the actual solutions. 

The main result is Theorem [7j which presents the two-scale variational problem for the 
functions h®(x), hf(x), h\(x,y), and hj(x,y) in the two-scale limits 

W{x) h°(x,y) = Xc(y)h° e (x) + Xi(y)h-(x), 

(15) 

VW{x) Xc(y)[V^(x) + V v hl(x,y)] + X i(y)[Vh?(x) + V y hl(x,y)], 
and its equivalence to the homogenized Maxwell system together with the unit cell problem 



(section 4.2) that determines the gradients of the corrector functions hi and h\ as well as 
the effective coefficients fi* and e*. 

In Theorem [8j we obtain the strong two-scale convergence of h v and Vh 71 to their two- 
scale limits. Strong two-scale convergence strengthens two-scale convergence by asserting 
convergence of energies. Theorem [9] asserts the convergence of the transmission and reflection 
by the micro-structured slab to that of the homogenized one. 

Acknowledgments 

The authors would like to acknowledge the inspiration of the IMA "Hot Topics" Workshop 
on Negative Index Materials in October of 2006 at the University of Minnesota, cosponsored 
by the Air Force Office of Scientific Research. This Workshop provided the foundation and 
stimulation for our work in this area. 

The effort of R.V.K. was supported by NSF grants DMS-0313890 and DMS-0313744, and 
the effort of S.P.S was supported by NSF grant DMS-0313890 during his visit at the Courant 
Institute and by NSF grant DMS-0505833. 



2 Overview of results 

This section serves to highlight the main results of the calculations of Section |4| including 
the "unit-cell problem" for the corrector functions and definitions of the effective electric 
and magnetic coefficients. In the presentation, we try to illuminate the physical meaning of 
the results and the methods of averaging. 

The feature of our system that leads to interesting magnetic behavior, as previously 
discussed, is the high surface conductivity of the resonators, which scales inversely to the 
length of the period micro-cell, that is, o~ 8 (x) = o~(x,x/r))/r) for x on the surfaces of the 
resonators. The dielectric and magnetic coefficients in the matrix exterior and interior to 
the surface of the resonators remain bounded, and are given by e{x,x/rj) and n(x,x/r)). 

The Maxwell system for if-polarized fields reduces to an elliptic equation for the scalar 
magnetic field h which includes interaction with a current on the surface of a conducting 
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resonator. Because of the extreme value of the conductivity, to leading order, the mag- 
netic field h behaves on the microscopic scale as a piecewise constant function, with jump 
discontinuities at the micro-resonator interfaces given by the surface current. The electric 
field, by contrast, exhibits fine-scale oscillations of order 1 in the matrix, and its tangential 
component is continuous across the micro-resonator interfaces. 

The exterior and interior values of h, as functions of the macroscopic variable x, are 
related to each other through a fixed complex ratio m(x) determined by the shape and con- 
ductivity of the micro-resonator, the magnetic coefficient in its interior, and the frequency. 
Representing the leading-order behavior of the exterior and interior values of h by the func- 



tions h° c (x) and h®(x), we will show in section 4.2 that 



h? 1 (x)=m(x)h° e (x), m(x)= -, / (x) -, , , (Hi) 

p[x) — iuofi[x) 

in which fx and p are weighted averages of the magnetic coefficient and the resistance (the 
reciprocal of the conductivity) of the micro-resonator. Thus as a function of macro- and 
microscopic variables, the leading behavior of h is given by 

h°(x,y) = M(x,y)h° e (x), (17) 

where M(x,y), which characterizes the microscopic variation of h, takes on two values: 
M(x,y) = Xg(u) + m ( x )XG* (y)- As we have mentioned in the introduction, the physical 
principle behind expression for m(x) is a balance of electromotive forces, which, in the 
microscopic variable, takes the form 

(h?(x) - h° e (x)) My) = M*, y)h-M <^(y) • (is) 



8G°( x ,y) JG 

The question now arises as to which value of h should be recognized as the appropriate 
macroscopic magnetic field in the homogenized bulk medium. The model problem of scat- 
tering by a slab that we have chosen represents a practical class of problems in which the 
boundary of the composite cuts through the dielectric matrix only so that the boundary con- 
ditions involve only the exterior H field. In our case, we have interface conditions between 
the slab and the surrounding air (continuity of h and E-t). Therefore we arrive at effective 
equations that govern h®(x) and the usual micro-cell average E av (x) of the E field: 

V L -E av {x)-iu^{x)hl{x) = 0, . 

(homogenized system) (19) 
V ± /i°(x) -iue*{x)E av {x) = 0. V ; V ; 

This system is evidently of two-dimensional Maxwell form, which may be seen a postiori 
as another justification for choosing the exterior value of h. Indeed, if the micro-resonators 



vary macroscopically, implying that m(x) is not constant, then the system (19), rewritten 
in terms of h® or some combination of and h®, say h® v (x) = a(x)h®(x), would contain an 
extra term — a(x) _1 V ± a(x)/i° v (x) in the second equation of ( JT9| , placing the system outside 
of the usual Maxwell type. The reason that this term is not present with the choice /i° lies 
in the unit cell problem for the corrector functions. 
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Because to leading order E-t vanishes on the surfaces of the micro-resonators, the cell 
problem is decoupled into parts exterior and interior to the domain G. It gives the order-?] 
corrector h l (x, y) as well as the leading-order part of the E field, E°(x, y). In its weak form, 
it reads 

/ E° c (x,y)-V ± v(y)dA(y) = for all v G Hg(G*), 

Jg* (exterior cell problem) (20) 

E° e (x,y) = ^(x^y)- 1 {V x h° e (x) +V y L h 1 e (x,y)) , 

in the exterior, which is a periodic-Neumann boundary- value problem for h\{x,y) and 
I E?(x,y)-V ± v{y)dA(y) = for aU v G H 1 (G) , 

Jg (interior cell problem) (21) 

E9(x,y) = ^(x^)- 1 (V^x) + V^hl(x,y)) , 

in the interior, which is a Neumann problem for h\{x,y). The interior cell problem has a 
simple explicit solution: 

V^h\{x,y) = -V^{x), . 

y (interior) (22) 

Ei{x,y) = 0. 



The second equation of (19) relates to E av (x) = f E°(x,y)dA(y). It is obtained by 



averaging E°(x, y) over Q, which is equivalent to integrating E®(x, y) over G* (since = 0), 
and using the definition in the second equation in (20). Thus E av (x) is a linear function of 
V^Or): 

E, v (x) = - [ s(x,y)- 1 (V ± h^x) + Xhl(x,y))dA(x) = -s*(x)- 1 \/ ± h c (x). (23) 
IU J G * iu 

This expression defines e*(x) in the usual way as an effective tensor for a periodic medium 
with inclusions that are perfect conductors, inside of which the electric field vanishes. In 



Section |4.4[ we argue that e* is an effective dielectric permittivity in the sense that _D av = 
£*E a _ v , where the average electric displacement field _D av is defined in the appropriate way. 
The effective magnetic permeability n*(x) is defined by 

H*{x)= [ M{x,y)fx(x,y)dA{y), (24) 
J Q 

and arises upon integrating over Q the equation 

V x x -E\x,y) + ■E\x 1 y) - iu ^(x,y)h°(x,y) = 0, for y $ dG, (25) 

where E l is the order-?] corrector to E°. 

The noteworthy feature in this definition is that, even if \x is real (and, say, constant 
throughout Q), fi* will still be complex- valued. This is because the usual cell average is used 
in averaging the B field fih whereas the average of h is taken exterior to the micro-structure. 
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The physical correctness of this method of averaging was noticed by Pendry, et. al. In section 
4.4[ we give a discussion of the averages of the four fields E, h, D, and b and how they emerge 



naturally in the effective equations. 

Finally, we come to the specific problem of the scattering of plane waves by a slab of the 
composite medium. Of course, one expects that in the quasi-static limit, the reflection and 
transmission coefficients will approach those of a slab of a homogeneous medium possessing 
the effective coefficients e* and /i*. In fact, in practice one would like to be able to deduce 
these bulk properties from scattering measurements. This has been the subject of many 
works in the literature, for example those of O'Brien and Pendry [T3J E] and Smith, et. ai, 
[22] . in which the scattering data obtained from a transfer matrix method for a periodic slab 
are compared to those obtained from expressions for the effective tensors of the corresponding 
homogeneous medium. In the present work, this agreement is established in the quasi-static 
limit by demonstrating that the limiting values of the electric and magnetic fields to the left 
and right of the slab are precisely those associated with scattering by a homogeneous slab. 
This result is embodied in the statement that the limiting electromagnetic fields outside the 
slab connect continuously to the fields E av and within the slab, where they solve the 
effective equations (19). 



3 Scattering by a slab with micro-resonators 

As our model problem, we choose the scattering of plane waves by a micro-structured slab of a 
composite material formed by periodically embedding metal micro-resonators in a dielectric 
matrix. The width of the slab is fixed, as the size of a cell tends to zero. We allow the 
material properties as well as the shape of the micro-resonators to vary on the macroscopic 
scale, but require that this variation be periodic with period 2ir in the x 2 direction, parallel 
to the slab. The fields will then be allowed to be pseudo-periodic in x 2 . 

As we have discussed, an important feature of the model is the requirement that, at the 
slab-air interface, only the dielectric matrix have contact with the air. 



3.1 Reduced Maxwell equations 

Due to the macroscopic periodicity of the structure and the pseudo-periodicity of the fields, 
the analysis of the scattering problem can be restricted to the strip 

S ■= {(an, x 2 ) ■ < x 2 < 2tt} . (26) 

As illustrated in Figure [TJ a period of the slab occupies a fixed subdomain Q of S, bounded 
below and above by the lines x 2 = and x 2 = 2tt and on the sides by XL (at x\ = a) and 
E + (at Xi = b). The union of the regions enclosed by the micro-resonators contained in 
Qo is called the interior domain and its boundary dfli is identified with the resonators 
themselves, or, perhaps more correctly, the surfaces of the resonators. 

Because we consider electromagnetic fields in a linear medium, we may treat harmonic 
fields at fixed frequency u>, 

£(x 1 ,x 2 ,x 3 ,t) = E(x 1 ,x 2 ,x 3 )e~ luJt , H{x 1 ,x 2 ,x 3 ,t) = ii(x 1 ,x 2 ,x 3 )e~ luJt , (27) 



12 



which, when inserted into the time-dependent Maxwell system, yield the harmonic Maxwell 
system for the spatial envelopes of the fields, 



VxE - icuuH = 0, 

28 

VxH + iueE- J = 0. v 1 

In this work, we study two-dimensional structures, in which the medium and the fields are 
invariant in one spatial dimension. In this situation, the Maxwell equations decouple into 
i?-polarized and iJ-polarized fields; we study the latter: 

H(x 1 ,x 2 ,x 3 ) = (0,0,h(x 1 ,x 2 )), 

E(x 1 ,x 2 ,x 3 ) = (Ei(x 1 ,x 2 ),E 2 ( y x 1 ,x 2 ),0), E(x 1 ,x 2 ) = (E\{xi, x 2 ), E 2 (x!, x 2 )), (29) 
J(x 1 ,x 2 ,x 3 ) = {Ji{x u x 2 ),J 2 {x u x 2 )^), J{x u x 2 ) = {Ji{x l ,x 2 ),J 2 {x 1 ,x 2 )). 

With the notation 

(the formal adjoint of V x - is — V x ), the curls of E and H are expressed as 

(31) 

vxn = —v~n, 
and the Maxwell system reduces to 

7- 1 . E - \ui u.h = 0. 

(32) 

The relevant fundamental theorem involving the operator V -1 - is 



VxE = V^-E, 
VxH = -V ± h. 



\7 ± -E-iufih = 0, 
V ± h-iueE + J = 0. 



V ± -FdA= / F-tds, (33) 

R JdR 

in which F is a vector field, R is a bounded domain in M 2 with boundary dR, and t is the 
counter-clockwise unit tangent vector to dR. 

We incorporate the regular part of the current J into the electric displacement field D = 
eE by allowing e to be complex)^} and retain only the singular part of J in the equation. As 



explained in section 1.3, our resonators are characterized by a complex surface conductivity 



a s relating the surface current to the electric field on the boundary of each resonator: 

J = a s (E-t)5 dUi t. (34) 



3 With e = e' + is", the conductivity in the material away from the metal rings is a :— cue" . In — \ujeE = 
—iioe'E + aE, the term —iuie'E is identified with the time derivative of the electric displacement and aE 
with the current density J off the surface of the resonator. 
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Fixing a Bloch wave vector k in the x 2 direction in the first Brillouin zone, k G [—1/2, 1/2), 
the scattered magnetic field h sc (total field minus incident field) has the outgoing form 

oo 

h sc (x 1 ,x 2 )= J2 c±e i{{m+K)X2+uM \ xi < a or n > 6, (35) 

m=— oo 

in which the exponents u m , which depend on to and defined by 

+ (m + k) 2 - Sq^quj 2 = (36) 

and the convention that v m > if > and iu m < if < 0. The incident magnetic 
field is 

h- mc (x 1 ,x 2 ) = e^+^e™, ^ 

in which > 0. 



3.2 Microstructure 

The micro-structure at each value of x is described by means of a microscopic variable 
y G 1R 2 and its fundamental period cube Q = [0, l] 2 (Fig. [2j. The boundaries of the micro- 
resonators are defined through a microscopic domain G in Q with boundary dG and exterior 
domain G* = Q\G and three complex- valued functions of both macroscopic and microscopic 
variables that have period Q in y. The domain G is assumed to be simply connected and 
contained wholly within the unit cell Q. Our analysis could be extended to the case in which 
G consists of more than one simply connected component or in which the micro-resonators 
cannot be modeled by a domain contained wholly within the unit cell. However, our model 
seems to include all of the split-ring resonators we have seen in the literature. 

The dielectric permeability e(x, y) is a tensor, whereas the magnetic permeability /i(x, y) 
and the conductivity a(x,y) are scalars. The tensor e, considered as a matrix, is symmetric, 
the real parts of all three quantities are positive and bounded from above and below, and 
their imaginary parts are semidefinite: 

0<e_ <Ree{x,y)£-Z<e + , < lme{x,y)£-£, x G Q , V e Q, f G R 2 , 

< IX- < Refj,(x,y) < < lmfx(x, y), x G f2 , y G Q, (38) 

< er_ < Rea(x,y) < a + , lma(x,y) < 0, x G fi , y G dG. 

Both e and \x are assumed to be continuous in x, and continuous in y off the boundary dG. 
Thus we allow different values in the interior and exterior of the resonators. To define the 
actual micro-structure at a fixed fine scale rj, we set 

e(x, x/rf), x G O , 

fj,(x,x/r}), x G Oq, (39) 
a(x,x/i]), x G dfli. 



fjP{x) = 
a^ix) = 
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As we wish to allow the conductivity of the metal cylinders to tend to infinity with 1/rj, we 
take 

a s = -a\x) (40) 
V 



as the surface conductivity in equation (34). Outside the slab, we take e and \i to be constant: 

e v (x) = e Q and fi n (x) = fi for x G Q\fi - (41) 

The interior domain Oi depends on 77 and is expressed in terms of the domains G through 

fl^fiefij: x/r] G G + n for some n G Z 2 }. (42) 

We have discussed the restriction that the edges of the slab not intersect the resonators. 
It is convenient, however, to assume a bit more: that the width of the slab encompass an 
integer number of period cells. Therefore we require rj = (b — a)/n for some integer n. Since 
we also assume that the structure is 27r-periodic in the slow variable x%, we also require that 
rj = 27r/m for some integer m. These conditions are equivalent to the condition that b — a 
is rationally related to tt. The set of permissible values of rj is denoted by T: 

rj e T. (43) 

For each rj, we let E v , h v denote a solution to the scattering problem, which is posed in 
its strong PDE form as follows. The subscripts e and % refer to exterior and interior values 
of functions. 

Problem 1 (Scattering by a slab, strong form) Find a function h v and a vector field 
E 71 in the strip S such that 



on 



S \ dni (44) 



V ± --B r? -iuj^hP = 
V ± h^ - lue^E* = 
E v 4 continuous across (45) 



hP e -hV + — E v -t = ondQ?, (46) 
W ( Xl , 2?r) = e iKX2 hP {xx , 0) and d X2 h r < ( Xl , 2vr) = e iKX2 d X2 h v {x u 0), (47) 

CO 

hP{x) = h inc {x 1 ,x 2 ) + a m e l((m+K)X2 ~ VmXl) for x x < a, (48) 

m=— oo 

oo 

W{x)= b m e^ m+R)x2+VmXl) for x x > b. (49) 

m=— oo 

Equivalently, one could pose the PDE as a divergence-form elliptic operator on h v alone and 
use the second Maxwell equation as the definition of E v : 

v-.(( £ .o-'v^ ) + . V /^o| 

EP = i (£")- 1 V ± /i" J 

The first equation reduces to V • {{e ri )^ 1 V h v ) + u 2 fi v h v = if e is scalar, and conditions 
involving E v -t are parsed in terms of h v through E v -t = A(fcxe" 1 V 1 /i I, )-n . 
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4 Formal expansion analysis 

In the formal asymptotic analysis of the strong form system, we assume that hP and E' n have 
expansions 

hP{x) = h°(x, x/rj) + rjh}{x, x/rj) + 0(r] 2 ), (51) 
E v (x) = E°(x, x/rj) + r]E\x, x/rj) + C(r/ 2 ), (52) 

and that these expansions can be differentiated term by term. By inserting these into the 
scattering Problem [TJ the results announced in Section [2] emerge. 

4.1 Expansion of the Maxwell system 

The first of the Maxwell equations, V ± -E v — iu fi v h v = 0, gives 

^V^-E^x.y) + (V^-E (x,y) + V/ -E^x, y)) - iufz(x, y)h°(x, y) = 0( V ), 

which yields the equations 

V y L -E°(x,y)=0, 

-E°(x, y) + Vy 1 -E\x, y) - iu fi(x, y)h°(x, y) = 0, 

Similarly, the second of the Maxwell equations, V~ L /i T? — iu eE = 0, gives 

r]~ l ^yh\x, y ) + (V£h°(x,y) + V^h l (x,y)) - iue{x,y)E°{x,y) = 0(r,), 

which yields the equations 

V y ± h°(x,y)=0, 

h°(x, y) + V/ h\x, y) - iu e(x, y)E°(x, y) = 0, 

The interface conditions tell us that 

E°-t and E 1 -t are continuous on dG, 
E°-t = 0, 

h° e -h? + aE 1 -t = 0, 

Let a micro-cell (Q + m)/r] with m G Z 2 depending on rj be chosen such that it contains a 
fixed point x, and set x = mrj. If we integrate in x over the part of Q{ and its boundary 
contained in this micro-cell, we obtain, after making the change of variable x' = x + rjy, 

r] [ ^{h n e - h?)ds{y) = -r) [ E v ■ t ds(y) = —rf f (V -1 -^) dA(y) = —iuif f fih v dA(y), 

JdG ° JdG JG JG 

( (55) 

in which the functions of x' and y are evaluated at (x + rjy, y). The expansions of E 11 and 
h v then yield the equations 

/ (j-\x,y)(h:(x,y) -h?(x,y))ds(y) = -iu; [ n{x,y)h?(x,y)dA{y) (56) 

JdG JG 

at each order n = 0, 1, 2, etc. 



for y <£_ dG. 



(53) 



for y ^ dG. 



(54) 



on dG. 
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4.2 The cell problem and the homogenized system 



From the first of the pair (54), V^~h°(x, y) = 0, which is valid for each i6ll and for y (jL dG, 
we infer exterior and interior values of the magnetic field that are independent of y: 



h°(x,y) 



'h° e (x), y e G*, 
h$(x), yeG. 



(57) 



The relation between the exterior and interior values is obtained from (56), which expresses 
a balance of electromotive forces on the boundary of a single inclusion, 



p(x)(h° e (x) - h$(x)) + xujjl{x)h^x) = 0, 



(5? 



in which 

M=[**,vWv), «*)=/ «*,vr*<v), (59) 

JG JdG 

and we obtain hi(x) = m(x)h e (x) and thus an expression for the magnetic field in a cell in 
terms of its value exterior to the inclusion, 



h°{x,y) = M(x,y)h° e (x), 



in which 



M(x,y) 



1, yeG*, 
m(x), yeG, 



m(x) 



p{x) 



p(x) — iujp(x) 



(60) 



(61) 



We next use the first equation of (53) and the second of (54), together with the continuity 



of E-t on the interfaces to obtain the unit cell problem 
V y ± -E°(x,y)=0, 

V^h°{x,y) + V^h 1 (x,y) -iuje(x,y)E°(x,y) = off dG, 
E Q e (x, y)-t = E?(x, y) ■ t = on dG, 
h x (x,y) and E (x,y) periodic in y. 



(62) 
(63) 
(64) 
(65) 



This problem determines h 1 and E° as functions of y, for each value of x. It is an inho- 
mogeneous periodic-Neumann problem with input V x h°, in which the exterior and interior 
first-order corrections, h\ and h\ are decoupled due to the homogeneous boundary condi- 
tion on the interface dG and are determined up to two additive constants in y, which are 



functions of x. One relation between these functions of x is provided again by (56): 



°" \x,y)(hl(x,y) - h}(x,y))ds(y) = -iu / p{x,y)h\{x,y)dA{y) 



dG 



(66) 



G 



As we have mentioned in Section [2j the interior value of E° vanishes identically. 
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To obtain a PDE governing the average exterior H field h® and the average E field E. 



we integrate the second equation of the pair (53) over the unit cell Q, as well as the second 



equation of (54) after applying e(x,y) 1 . The result is the homogenized system 

V-L-EnW-hjfSWhKx) = 0, 

V^Or) - \ue*{x)E ay {x) = 0. ^ ' 

The effective magnetic permeability tensor /i* is given by 

H*(x)= [ M{x,y)fi{x,y)dA{y), (68) 
Jq 

and the effective dielectric tensor e* arises as we have explained in Section [2] 

P av (x) = / E°(x,y)dA(y) 
Jq 

= -[ e(x, y )- 1 {V ± h e (x) + V y ± h 1 c (x,y))dA(y) = -e*(x)- 1 V ± h° e (x). (69) 
vjj Jg* ioj 

In order to justify calling e* an effective dielectric permittivity, we must ensure that e*E 
is correctly interpreted as an appropriate cell average of the electric displacement field D. 
We address this issue in section 14.41 



4.3 Corrector functions 

Because of the condition (64) that E®(x,y) ■ t = Ef(x,y) ■ t = on dG, the interior and 
exterior gradients of the corrector function h 1 (x, y) are decoupled and are determined as 
linear functions of Vh®(x) and V/if(x) through the cell problem restricted to G* and G, 

^hl(x,y) = P e (x, 2 /)V ± /i e °(x) in G*, (70) 
V^h\{x,y) = P^y^h^x) in G. (71) 

The corrector matrices P e (resp. Pi) is defined by setting P e £ (resp. p£) equal to the unique 
solution of the cell problem in which £ stands for V h®(x) (resp. Vh®(x)). Posed in their 
weak forms, these problems are 

/ e(x,y)- 1 ^ + P c (x,y)^V ± v(y)dA(y) = for all v G H#(G*), (72) 
Jg* 

[ e(x,y)- 1 lZ + P i (x,y)Z]-V ± v(y)dA(y)=0 for all v G H l {G). (73) 
Jg 

As we have seen, p admits a very simple form: 

Pi(x, = (74) 
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The tensor e*(x) is defined through 



Jg* 



(75) 



The term of inhomogeneity in the cell problem, v*^~ h°(x, y), might as well be parsed in 
terms of the average exterior H field h° e and its gradient, 

* V) - V- (M(x, = {™ + v±m{m (76) 

From this point of view, the corrector function h l {x,y) is a linear function of V h®(x) and 
h®(x) for any fixed value of x, so that it is determined in two parts, given by a corrector matrix 
P(x, y) applied to the gradient V _L /ig(x) and a corrector vector Q(x, y), which multiplies the 
scalar h®(x): 

V^h\x,y) = ^.yjV^^+^i/)^), (77) 
and one computes that these are related to P e and Pi by 

P(x, y) = Xc(y)Pe{.x } y) - Xi(y)m(x), (78) 
Q{x,y) =Xi(y)V ± m(x). (79) 



4.4 Discussion of average fields 



The homogenized system (67) is clearly of Maxwell form, in which E av (x) and h®(x) play the 
role of electric and magnetic fields in a bulk two-dimensional medium. These fields are certain 
micro-cell averages of the actual electromagnetic fields in the micro-scale composite, E av (x) 
being the usual cell average and h° e (x) being the average over the part of the cell exterior 
to the micro-resonator. We expect additionally that the fields e*(x)E av (x) and fj,*(x)h®(x) 
should represent suitable cell averages of the D field eE and the B field (0,0, fj,h). 

In the medium exterior to the slab (air, for example) all of these fields reduce to E, h, 
D = eE, and b = [ih. In the composite, one has be careful about how cell averages are 
to be understood. The standard point of view in the physics literature (see, for example, 
[TT] or [19]) is based on the fact one should keep in mind the different geometrical roles of 
these fields: Averages of the E- and H fields should be taken over line segments traversing 
a micro-period in the direction of each field component separately because these fields are 
naturally integrated over one-dimensional curves (they are fundamentally one- forms). On 
the other hand, averages of the D- and B fields should be taken over sides of the unit cube 
perpendicular to each component separately because these fields are naturally integrated 
over surfaces (they are fundamentally two-forms). 

In the two-dimensional if-polarized reduction, this scheme amounts to computing aver- 
ages in the following way. Since H is out of plane, its average is taken simply by evaluating 
h at a suitable point, whereas the average of b must be taken over the entire unit cell. The 
Xi-component of the E field should be averaged over line segment in the x^-direction, whereas 
the average of the Xj-component of the D field reduces to the average over a line segment 
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in the ^/-direction, where i' — (i + l)(mod2), because D is constant in the out-of-plane 
direction. 

We now demonstrate that, from this point of view, the tensors e*(x) and fi*(x) that 
emerge in the homogenized equations do indeed relate the natural cell averages of E and D 
to each other as well as those of h and b. 

In defining /i av , a choice between its evaluation in the exterior or the interior of the 
resonator needs to be made. In problems of scattering by a slab, one avoids cutting through 
a micro-resonator and therefore arranges the slab-air interface such that the air is incident 
with the dielectric matrix of the composite that is exterior to the resonators. The exterior H 
field is therefore the one that connects continuously with the H field in the air, and therefore 
we should work with the average of h in the exterior: 

h av (x) := h° c (x). (80) 

Evidently, since h is constant in y to leading order exterior to the resonator and constant in 
the out-of-plane direction, no averaging of oscillations is truly taking place in this definition. 

The B field, which is identified with its out-of-plane component b = fih, should be averaged 
over the cell Q: 

bUx)-= f Kx,y)h°(x,y)dA(y), (81) 



which is exactly equal to (iu) 1 times the second term in the right-hand side of the first 



equation in (67), and we have therefore 

K w {x) = /j,*(x)h av (x). (82) 

Now let us examine the E and D fields. Let e\ and €2 denote the vectors (1, 0) and (0, 1), 
and set i' = (i + l)(mod2). The cell average of the e« component of the E field should be 
taken over a line segment in the direction of traversing a period cell, 

i 

E &v {x)-ei:= J E°{x,y)-e i dy i . (83) 

Vi = 
Vi> = c 

These integrals are independent of y^ because -E°(x,y) = 0, and therefore by this 
definition, E av can be taken to be the area integral of E°(x, y) over Q, as we have defined it 



above in (69). 



The 6i component of the D field D = eE is averaged over a line segment in the direction 



of e,-/ 



D av (x)-ei := / e(x,y)E°(x,y)-e i dy i/ . (84) 



w = o 

Hi = c 
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If the path of integration remains exterior to G : this yields 



i 



D av (x)-e l = - I (V a ± /i e °(x) + V^ 1 (x, 2 /))-e l d^ = -V :c ± /i e (x)-e i , (85) 
iuj J y iu; 



Vi = c 



and therefore, with this definition of D av , we obtain from the homogenized equation (67) 



D^{x) = e*(x)E av (x). (86) 

5 Two- scale convergence analysis 

This Section establishes with mathematical rigor the results of the formal analysis. Our 
main result is that the solution of the problem of scattering by the micro-structured slab 
tends to the solution of the problem of scattering by a homogeneous slab. This amounts 
to convergence of the electromagnetic fields to the average fields discussed in the previous 



section and the fact that these average fields satisfy the effective system (67) in the slab. 
Of course we must also show that the field /z av = h° e connects continuously to the function 
h in the air and that the tangential component of E av at the edge of the slab connects 
continuously to that of E from the side of the air. 

These things are properly handled by means of the weak form of the scattering problem. 
The outgoing conditions for the scattered field are most conveniently expressed by employing 
the Dirichlet-to-Neumann map T on the sides T± for outgoing fields. This means that fields 



of the form (35) are characterized equivalently by 

dh 



dn 



r± = —T (h sc \ r± ) . (87) 



i 

For g e Hz(T±), and {g^} denoting the Fourier transform of ge~ 1KX2 , this map is defined 
through 

{TgT m = . (88) 

i _i 

T is a bounded operator from H2(T±) to H 2 (F-|-), with norm \\T\\ H1/2 _> H -i/2 — C\ + 
C 2 |o;| for some positive constants C\ and C%. For real values of u and k, it has a finite- 
dimensional negative imaginary subspace spanned by the functions e l( - m+K ' ,X2 for which > 
and a positive subspace equal to the space orthogonal to the imaginary one. We denote the 
corresponding decomposition of T by T = T + + iT~. 

Denote by H^D^) the K-pseudoperiodic subspace of if 1 (fi^)- 



#M = e H^OH) : u\ X2=27T = e 2 ™u\ X2=0 } ■ (89) 

It is natural to identify H^iPQ) with a subspace of L 2 (Q) by extension by zero into and 
H 1 ^ 1 ?) with a subspace of L 2 (Q) by extension by zero into O^- We can then define the space 
V v in which h v resides by 

H 1 ^) © H\n?) =: V v ^ L 2 (Q) ^ (V 11 )* V n . (90) 
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It will be convenient also to define V° := H^(Q). 

The following weak form of the scattering problem is equivalent to the strong form for 
smooth s, /i, and a and for functions h e and hi that are smooth in Q e and Qi . 

Problem 2 (Scattering by a slab, weak form) Find a function h% = h v = h^(Bh^ G V v 
such that 

[ [((e")- 1 V 1 /i")-V I i;-wVH dA-iu I rj^y 1 ^ - h?)(v e - v$ ds + 
Jq L -I JdW 

+ e 1 J (Th2)v e dx 2 = -2i^£o 1 J e ^ m+K > 2+1/ ^v e dx 2 (91) 

for all v = v c © Vi G V v , and let E v G L 2 (Q) be defined by 

E v = — (e^V^hP in ft!? U fi?. (92) 
iu; 

This scattering problem always has a solution. It is possible to prove that on finite 
intervals of the real a;- axis avoiding a discrete set of frequencies, is unique for all 77 
sufficiently small, but this fact will not needed in our analysis. 

Theorem 3 (Existence of solutions) For each frequency uo and rj G T, the scattering 
Problem^ has a solution h v . 

Before proving the Theorem, we look more carefully at the forms involved in the weak 
formulation of the scattering problem. 

Let the form a^(u,v) in V v be defined by the left-hand side of the equality in Problem 
[2j with u in place of h n ' . We split into two parts: 

al(u,v) = bl(u,v)-aj 2 c^u,v), (93) 



bl(u,v) 



[ ((e v )- 1 V ± u)-V ± vdA + (94) 
Jn 



-iuj rj(a v ) 1 (u c - Ui)(v e - v{) ds + e 1 / (Tu e )v e dx 2 , 
Jan"? Jt ± 

c v {u,v) := [ yPuvAA. (95) 

The incident plane-wave field provides the forcing function for the scattering problem, given 
by the functional / on V v defined by 

f(v) := -2i^£- 1 I e^ ifl+ ^ X2+v ^v e dx 2 . (96) 



The form b^ is coercive, with constants that are independent of rj, and / is uniformly 
bounded: 

K(u,v)\ < ( Tl + (97) 

Reffl(u,u) > 5\\uf vr , , (98) 

\\f\\ (v , r <C. (99) 
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The second inequality holds because of the upper bound on e v . In order to establish the upper 
bound on \a v (u,v)\, one uses the bound on T discussed above and the following Lemma [4 
showing that the term involving the interface dfl^ is a bounded form on V v . The bound on 
/ is proved in Lemma [5] 

Lemma 4 There exists a constant C such that the following hold for each rj £ T. 
1. For all u = u c © tij and v — v e © v± in V v , the following inequality holds: 



an? 



T)(o£) (tie - Ui)(v e - Vi)ds 



< C\\u\ 



V \\ v II Vi ■ 



(100) 



2. For each u = u c © tii £ V 7 , t/iere are functions u e £ H^ft) and U{ £ if 1 (fi) suc/i £/ia£ 
w e In? = u e \ n n and «i| n ? = itiln? 



Iwelln 1 ^) — Clplljfi^)' ll ?i ill// i (n) - C|l n lli/i(n;')- 



(101) 



Proof. To prove part (1), we use a Poincare inequality in the square Q: For w c £ H (G*) 



My)| 2 d*(i0 <C (\w e (y)\ 2 +\Vw e (y)\ 2 )dA(y). 

8G JG* 



(102) 



Rescaling this estimate to a cell of size r\ in the variable x and summing over all cells in Q 
yields 



\u e (x)\ 2 ds(x) = — [ (|m c (x)| 2 + ?7 2 |VMe(a;)| 2 ) dA(x). 
any V Jn% 

This, together with an analogous estimate involving give 



I \ 2 A <- ^11 II 2 
tl e °- s _ — Pc mm' 



aw 



■ ■ e||^l(nj), J \Ui\ 2 ds < —\\Ui\\ 2 H1(Q vy 



(103) 



(104) 



Finally, recalling that cr_ < Reo~(x,y), we obtain 



an v 



rj(a^) (ti e - tii)(f e - Vi)ds 



\u e — Ui\ ds / \v e — t?i| ds 
anv Janv 



< C 2 (ct_) 2 [ l|ti e || H i (n g) + ||tii|| 

Hi(Q v ) I I II^bII mfn'n + llfi || m/r>i\ T/n t> 



\-2|L.I|2 I, |,2 



(105) 



An extension of ti = u e in part (2) can be obtained from the extension operator defined 
in [8]. For the convenience of the reader and coherence of this text, we include a proof. The 
extension is accomplished by defining u e to be a harmonic function in whose trace on dQ^ 
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coincides with that of u e . The extension is K-pseudoperiodic because u e is and dQ PI dfy = 0. 
Let x n = nrj G Qo where n G Z 2 . Using the following standard inequalities in the unit cube, 

\V y u c {x + V y)\ 2 dA{y) < c \V y u e (x + V y)\ 2 dA(y), (106) 

JO. JG* 

\u e (x + r]y)\ 2 dA(y) < c (\u c (x + r]y)\ 2 + \V v u c (x + r]y)\ 2 ) dA(y), (107) 

JG* 

we obtain the estimates in the scaled cubes x n + r\G* , 

/ \V x u c (x + x)\ 2 dA(x) < c / \V x u e (x + x)\ 2 dA(x) (108) 

Jr/Q JrjG* 



and 



/ \u c (x + x)\ 2 dA(x) <= c (\u e (x + x)\ 2 + 7] 2 \V x u e (x + x)\ 2 )dA(x). (109) 

JvQ JvG* 



>vQ JvG 
Summing up over all micro-cells in Q , we obtain 



/ \V x u e \ 2 dA(x) < c I \V x u e \ 2 dA(x), (110) 
I \u e \ 2 dA(x) < c I {\u c \ 2 + rj 2 \V x u c \ 2 )dA{x), (111) 



The result trivially extends to Q because the edges of the slab E± do not intersect the 
resonators. The extension of U\ is handled identically in fl and then extended to a function 
in H 1 ^) by a harmonic function in whose trace coincides with u~i on E±. The extension 
to Q can be taken to be in H^(D,) if ii\ G iy*(f2o)- ■ 



Lemma 5 Lei / G (V)* be defined as in (96). There exists a constant C such that 
\\f\\ {Vny <C for all TjeT. 

Proof. 

Let u = u e © Ui G V n be given, and let w e be the extension of w e provided by part (2) of 
Lemma HI Then 

1/0)1 = l/(OI = |/(w e )| < H/|lHi(n)«ll«ellHi(n) ^ c ll/lln-i(n)* IKIIhW) ^ c ll/llffi(n)*ll w IU' 

(112) 

which implies 

\\f\vny < c||/||Hi (n) .. (113) 

■ 

The forms we have defined give rise to bounded operators £*J and C" 7 from V v into itself 
such that 

bZ(u,v) = (B*u,v) V n , (114) 

C "(u,u) = (C%,u)y,, (115) 
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where (•, is the usual inner product in V v . Let / denote that element of V v such that 
(/; ')vn = /■ The scattering Problem M now takes the form 



al(u, v) = f{v) for all v G V v , or B*u - u^C^u = f. (116) 

We now prove that the scattering Problem [2] has a solution. We refer to Bonnet-Ben 
Dhia and Starling [3], Theorem 3.1, where this is shown for a similar problem of scattering 
by a slab. 

Proof of Theorem [3[ Since is coercive, is a bijection with bounded inverse, 
and since \x G L°°(0), C v is compact. Therefore a Fredholm alternative is applicable: The 
scattering problem ( |116 ) has a solution u if and only if (/, v) = for all v G Null(.B^ — ^C^Y , 
or equivalently 

/(f) = for all v E V v such that d^(w, v) — for all wj G V^. (solvability condition) 

(117) 

Now any function v satisfying the adjoint eigenvalue condition a^(w,v) = for all w G V v 
satisfies, in particular 

al(v,v) = 0, (118) 
and, since the imaginary part of a% is nonpositive, we find that 

r_(v) = o. (119) 

Therefore the trace of v on T_ possesses only decaying Fourier harmonics, exterior to the 
slab, 

v( x ) = 7 ±e i((m+K)x2+!yra|xi|) exterior to the slab, (120) 

m:iv m <0 

whereas / is defined through integration against the trace of the incident field, which pos- 
sesses only a propagating harmonic. Thus we obtain 

f{v) = -2iu M ! e i((m+«>2+^i) j ^ e l{ - {m+K)X2+VmXl) \ dx 2 = 0, (121) 



i.e., the solvability condition (117) is valid. Thus, by the Fredholm alternative, there exists 
a solution to 

al(hl, v) = f(v) for all v G V v . (122) 

■ 

Since we do not have a priori bounds on the solutions h v , we scale for the time being the 
magnetic and electric fields by a number mP G [0, 1], 

m" = min{l, \\h"\\j}}, (123) 
u T '(x) =m r 'h T i{x) 1 (124) 
F r >(x) = m v E v (x) } (125) 
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such that the vP are solutions of the scattering problem with a scaled incident wave and are 
bounded in the L 2 norm uniformly in rj: 

al(u v , v) = mPf{v) for all v G V\ (126) 

F v = — (^rVV 7 , (127) 

\\u v \\ L 2 < 1. (128) 

After proving the two-scale convergence of the u v to a solution of the homogenized equations 
(67), we will conclude that the hP are in fact uniformly bounded in V' n (at the end of the proof 



of Theorem |7j). This result uses the fact that the scattering problem for the homogeneous 
slab has a unique solution. It will be shown that this always holds true in the case that yP 
is real, in particular in the typical case of nonmagnetic materials, yP = y^. 

The next Lemma establishes the existence of two-scale limit^j of subsequences of vP and 
Vu v whose y-dependence reflects the discontinuities at the surfaces of the micro-resonators. 
No other feature of the weak-form Problem [2] is used here. The symbol "— ^" means "two- 
scale converges to" . We denote the characteristic function of G extended periodically to 1R 2 
by Xi(y) and set Xeiy) = 1 - Xi(y)- 

Lemma 6 (Two-scale limits) Every subsequence of T admits a subsequence T' and func- 
tions u° e (x) G Hl(n),uf(x) G H^Qq), ul(x,y) G L 2 (Q;H^(Q\D)/R), and u{(x,y) G 
L 2 (Q ;H^(Q)/R) such that, in Q, 

u"(x) Xe(y) u e( x ) + Xi(y)uf(x) ) 

> 71 G T' (129) 

W(x) X e(y)[Vu° e {x) + V y ul(x,y)} +Xi(y)[Vu°(x) + V y u\{x,y)]\ ' 



Proof. With v — vP in (126), we obtain 

b(vP, vP) = ivPf{vP) + u 2 c{u\ vP). (130) 

Then, using the coerciveness of b%, the boundedness of y, and part (b) of Lemma (J4J), we 
obtain the estimate 

S\\u% v < Reb(uP,vP) < \ f{vP)\ +u 2 [ Re//V| 2 dA(x) 

Jn 

< C\\u\\ vn + u 2 y + \\u v \\ 2 L2 <(C + u 2 y + )\\u ri \\ vv , (131) 

from which it follows that 

Kllv„ < (C + u 2 y + )/6. (132) 

We now have that u v and Vu v (defined for points off dQ?) are bounded sequences in 
L 2 (Q), and therefore, by Theorem 1.2 of Allaire p], there exist functions u° G L 2 (Q x Q) 
and £° 6 L 2 (fl x Q) 2 and a subsequence T'cT such that 

vP{x) u°(x,y), and W(x) £ (x,y), r] G T'. (133) 

4 A sequence vP{x) two-scale converges to u°(x,y) (u v u°) in if lim J( ^o Jq u v (x)(j>(x, x/rj) dA(x) — 
In Iq u °i x i y)4>{ x i y) d^4(y) dA(x) for all smooth functions <f>(x, y) in f2 x K 2 , periodic in y. 
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The behavior of these functions is well known to be trivial in the region Q \ Q outside the 
slab, where they have no y-dependence, and u°, which we denote by u® in this region, is of 
class H 1 there. 

Let * G Cg°(fio; C#(R 2 )) 2 with ty(x,y)-n = for y G dG be given. Upon integration by 
parts, the latter condition eliminates integrals over dfl? which would otherwise arise due to 
the discontinuity in vP on the interfaces of the resonators: 

J Vu n (x)-^(x,x/r])dx + J vP{x) (v x -^{x,x/r]) + ^V y -^(x, x/rj)j dx = 0. (134) 

Multiplying this equality by rj and using the two-scale convergence of u v and Vu v , the limit 
for r] G T' gives 

/ u°(x,y)V y -V(x,y)dydx = 0. (135) 
Jn Jq 

Let ip G C£°(Qo) be given and use a separable function (f>(x)$(y), with &(y)-n = on dG, 



in place of $?(x, y) in (135) to obtain 



4>(x) / u°(x,y)V-$(y) dydx = 0, 
Jo 



(136) 



from which we obtain 



u°(x,y)Vy$(y) dy = 0, 



(137) 



for each such It follows from this and the connectedness of G and G*, that u°(x,y) is 

independent of y in each of these domains. Therefore, for x G Qq, 



u°(x,y) = Xc(y)u° c (x) + Xi(y)u-(x) 
u2(x) Xe(2/)«°(a0, 
u?(x) Xi(2/)w°(a;). 



(138) 
(139) 
(140) 



To prove that u®(x) G H 1 ^^), we must prove that there is a constant such that, for each 
smooth function Q(x) with compact support in f2 , 



u° e (x)V-Q(x)dx 



9-o 



< const. 1101 



:i4ii 



By Lemma 2.10 of pQ, we can further impose upon \l/ (in addition to $?(x,y)-n = on dG) 

2 (Qq 



the three properties (i) V y -^(x,y) = 0, (ii) J G , *&(x,y)dy = Q(x), and (iii) l^Hia^oxG*) 2 — 



c||6|| L 2(q )2- One one hand, we obtain 



V vJl{x) {x , x / rf)dx = I u2(x)'V x -^f(x 1 x/r])dx — > / / u^(x)V x -^(x,y)dydx 

Sao J J £7o J 

= I u° c {x)V-[ ty{x,y)dydx= [ u° e (x) V-0(x)dx, r\ G T'. (142) 

Jfln Jg* Jq 
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On the other hand, 



/ Vu v c (x)-^/(x,x/r])dx — > / / £?(x,y)-^(x,y)dydx, r] G V , (143) 
Jn Jn Jg* 



so that 



fin 



u° e (x)V -<d(x)dx = / / ^°(x,y)-^(x 1 y)dydx, r] G T', 

Jfio JG* 



(144) 



and the inequality 



Co JG 



£°(x,y)-^(x,y)dydx 



< \\e\ 



< 4e\ 



101 



L 2 (n xG*) 2 ll w llL 2 (O ) 2 

(145) 

demonstrates that acts as a bounded linear functional on L 2 (Qq) 2 so that u° e G H 1 {VLq). 
An analogous argument proves that uf G H 1 ^^). One can see that itg and are in HI(Qq) 
by applying the arguments above to an extension of Qq into a neighborhood of X2 = 2n with 
u'' extended pseudo-periodically into this region. 



To obtain the functions ul(x, y) and u}(x, y) in Q , we use lim^o J a Vu^ix}^ (x , x / r])d 



x 



Jq >/g* ^^eC^) ' ^(x,y)dydx to obtain J* nn J G » (£°(x,?/) — Vm°(x)) -\I/(x, y)dydx = 0, which 



implies, as before, that 



m j g 



G* 



(^(x,y) - Vul(x))-$(y)dydx = 



(146) 



for all $(y) G H^(Q) 2 with V-$ = and $-n = on <9G, from which we infer the existence 
of a function u\(x,y) G L 2 (f2 ; H#(G*)) such that 

e°(x,y) = V Me °(x) + V^(a;, 2 /), y E G*. (147) 

In an analogous manner, one establishes the existence of u}(x, y) G L 2 (Q Q ; H^(G)) such that 

Z°(x,y)=Vu°(x) + V y ut(x,y), y G G. (148) 

To prove that u° e G i/ 1 (fi), we must prove that the trace of from the right and left of the 
slab boundaries £± are equal. By Lemma [4], the functions ul can be extended to functions 



vP G H 1 ^) in such a way that 



u' 1 



< C. We may then extract a subsequence that 



weakly converges in if 1 (fi), and by Proposition 1.14 (i) of [1J, this subsequence two-scale 
converges to its weak limit, which must evidently be equal to u®. ■ 

The key to obtaining two-scale convergence of the h v to a solution of the scattering 
problem for a homogeneous slab is the uniqueness of this solution. The weak form of this 



problem is given by (153). If Im(/i*(x)) > 0, then uniqueness is guaranteed. As noted in 
section 1.3, equation (13), Im(/i*(x)) is indeed positive provided that fi v (x,y) is pointwise 
real and positive. TheoremsVA Uft anduA require as implicit assumptions that Im(jj*(x)) > 0, 



or, more generally, that (153) admits a unique solution 
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Theorem 7 (Variational problem for the two-scale limit) The two-scale limiting func- 
tions h®(x), hf(x), hl(x,y), and h\{x,y) satisfy the following variational problem: 



nJG* 



(e(x, y)- 1 [V ± / i °(x) + hl(x, y)}) ■ [V Wfr) + V^v c \x, y)]dA(y)dA(x)+ 



[ [ (efayy^hKx) + V^!(x,|/)])-[VV(x) + V y ± v i 1 (x,y)]dA(y)dA(x)+ 
Jn Jg 

-uo 2 f If pi(x,y)dA(y)\ h° e (x)v e °(x)dA(x) - uo 2 f If pt(x : y)dA(y) 
Jn Ug* J Jn Ug 



h°(x)v i °(x)dA(x)+ 



lu / a(x,y) Ms(t/) 
Jn UdG 



(h° e (x) - h^(x))(v e (x)-v i °(x))dA(x) 

= -2iv fh £o 1 J e^+^+vmxi)^^ (149) 



for all v° c (x), vf(x) G H^Q) and v l c {x,y) G L 2 (fi ; H^(Q\D)) © L 2 (n \ Q ), v\{x,y) G 
L 2 (Sl Q -Hl(G)). 

This problem is equivalent to the the system 



h®(x) = m(x)h®(x), 

V£h\x,y) = xMPc{x,v)V L tii{x) + Xi{y)Pi{x,y)V L h!l{x) 
E°(x,y) = j^efoy)" 1 [Vh°(x) + V y ± h 1 (x,y)} , 

E^{x) = -e^xy^h^x), 
iuj 



[iuE^(x)-VH{x) - u 2 ^{x)til(x)v(x)] dA(x) 



(150) 
(151) 

(152) 



(153) 



= -2wrn£o 1 J e^ fh+K)x2+m) vdx 2 for all v G if 1 (ft), 
in which the average E field E w is defined as 

E av (x) := f E°(x,y)dA(y). (154) 
Jq 

Proof. Let T' be a subsequence of T such that mP converges, say to m° G [0, 1]. We use 
test functions of the form 



v ( x ) =Xc(x/v) [ v c ( x ) + V v l( x > x / r l)] +Xi( x /v) [vi(x) + r]vl(x,x/r])~\ = v°(x) + r qv 1 (x, x/rj). 

(155) 
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that are smooth off the boundary of the resonators in the weak-form Problem [2j 



Jn 

— to 2 I fi(x,x/rj)u v (x) (v°(x) + rjv x {x, x/rj)} dA + 
Jn 

T]cr~ (x, x/rj) (u2(x) — u^(x)) (v e (x) — Vi(x) + r](v e (x, x/rj) — Vi (x, x/rj))\ ds = 

= -2iv m J e l{{m+K)x2+v ^ Xl) (v e °(x) + r]v e \x,x/ri)) dx 2 . (156) 

The following analysis uses two-scale convergence results of pQ; we refer in particular to 
the proof of Theorem 2.3 of that work. Consider first the first term of (156). By assumption, 
Xe,i(y)^(x, y)^ 1 G C[f2 ; L^(IR) 2 ] 4 , and therefore we also have 



( e (x,y) := xM^fav)' 1 )* (VVW + V^-e 1 ^,^)) e C[n ;L c #K ,~,j s , 
d(x,y) := X i(y){e{x,y)~y (W(i) +^v^(x,x/r l )) G C[Q„;I|(Iff, (157) 
fey) = (e(x,i/)- 1 )t (V^ c °(x) + Vj-vjfax/ri)) G C[Q\Q ; L#(R) 2 ] 2 . 

This implies that the fields £ e ,i and £ strongly two-scale converge in their respective domains]^] 
As V J *w' ? two-scale converges in Qq, Theorem 1.8 of [lj^jjustifies passing to the two-scale limit 
in (156) in the entire domain Q. A similar argument applies to the second term in ( 156[ ). 

One deals with the limit of the interface term by applying the theory of two-scale conver- 
gence on periodic interfaces. We refer the reader to [2] and p2] for results on this subject. 
Proposition 2.6 of j2] applied to the extensions of and ul from Lemma [4] guarantee the 
two-scale convergence of the traces of these functions to u®(x) and u®(x) If a°(x, y) is taken 
to be continuous, then by extending it to a continuous function on x dG and incorporating 
it into the test function, we see that o" _1 (x, x/rj)(u2(x) — u?(x)), defined on dfl?, two-scale 
converges to a™ 1 (x, y)(u® (x) — u®(x)). 

The limit for rj G T' yields the two-scale weak-form PDE system for the two-scale limits 



5 Strong two-scale convergence of w v (x) to w°(x,y) means that w v two-scale converges to to and that 
/ \w v (x)\ 2 dx converges to // \w°(x, y)\ 2 dxdy. 

6 Essentially, strong two-scale convergence of u n (x) to u°(x,y) and (weak) two-scale convergence of v v (x) 
to v°(x, y) implies weak convergence of u v (x)v v (x) to J Q u°(x, y)v°(x, y)dA(y). 

7 A sequence w v (x) £ L 2 (dti?) with y J gQr , \w n (x)\ 2 ds(x) < C two-scale converges to w°(x,y) £ L 2 (Q x 
dG) if lim^o y J gn ^ w' n (x)(f)(x, x/rj) ds(x) — f n f gG w(x,y)(f)(x,y)ds(y)dA(x) for all continuous 4>{x, y) 
defined on Cl x dG. 
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tig, U®, Mg, and u\: 



[ [ [e(x,yr 1 (V ± u c (x) + V y ± ul(x,y))-(V ± v c (x) + ^v c 1 (x,y))}dydx + 
Jn Jg* 

+ [ [ [e(x,yr 1 (V ± u^x) + V y ± ul(x,y))-(V ± v i °(x) + V y ± v i \x,y))}dydx + 
Jn Jg 

— lu 2 u®(x)v e °(x) / jjt{x,y)dy dx — uj 2 / u®(x)vi°(x) / [i(x,y)dydx + 



G" 



G 



iuo / — M?(a;)) (f e °(x) — / a 1 (x 1 y)dydx 

Jn JdG 



= -2lVrn J e l{{fh+K)X2+m) V e °{x)dx 2) 

for all v e °, Ui° e H\n), v] e L 2 {Q- H^{G*)), v\ e L 2 {Q- H l # (G)). (158) 

Next we pull the two-scale variational problem apart to reveal the equations for the 
average fields and the cell problem for the corrector function. In Q , we first put v = and 
obtain 

/ / [e{x > y)-^V ± ul{x)+V^ul(x i y))-V^v e 1 {x > y)]dydx + 
Jn Jg* 

[ [ [e{x,y)~ 1 (V ± u9{x)+V^ul{x,y))-V y L v i 1 {x,y)]dydx = 
Jn Jg 

for all vl E L 2 {n- H^(G*)), v\ E L 2 (ft; H#(G)). (159) 



+ 



By using separable test functions v\{x,y) i— > <p(x)vl(y) and v\{x,y) \— > <p(x)vl(y), we obtain 
the cell problem for each x: 



L 



[eix.yY 1 (V x u° e {x) + V£u\{x,y)) ■ ^v e \y)] dy + 

+ f [eix^y 1 (VV(*) + V^ul&y)) ■ V^fo)] dy = 

for all vl e vl G H l # (G). (160) 



G 



Using this together with the definitions (72-73) of the corrector matrices P e and Pi, we 



obtain equation ( 151 ). In addition to v = 0, let us also set v I = and v\{y) = £-y to obtain 



/ e(x, y)- 1 (V ± w?(x) + u?(x, y)) ■ k x f dy 
Jg 



(161) 



It follows from this and (152) (taken as the definition of E°(x, y)) that 



E°{x,y)dA{y) = 



(162) 



G 
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Now let v®(y) be arbitrary and set v® = and v = 0. The gradient term vanishes by (161) 
and we obtain 



p{x,y)dy 



UG 



iq (x)v{ (x) dx + 



+ iu 



a(x,y) l dy 



n UdG 



or 



u?(x) 



p(x) 



(u° e (x) - u?(x)) iTi°(x) dx = 0, (163) 

(164) 



f)(x) — \LOp{x) 

which, by definition (61) of m, is equation ( |150 ). Now set v\ = and let vf and v® be 
arbitrary to obtain 



/ e{x,y)- 1 (Vu°{x) + \/^u 1 (x,y))dy 
Jo. 



V ± v°(x)dx + 



— UJ 



p{x,y)M(x,y)dy 



u° e (x)v°(x)dx = f(v°). (165) 



This, together with the definition (154) of E av (x) and equation (162) gives 



E°{x,y)dA{y), 



(166) 



G* 



and using the definition (75) of e*(x) and the definition (68) of p*(x), we obtain (153). We 



observe that equations (152-153) constitute the weak form of the homogenized system (67). 



the two-scale weak-form system (158). 



One also confirms that equations (150 -151) together with the definitions of e* and p* imply 



Now, we must prove that ||/i' ? || i 2 < C for all r\ G T' so that u v can be replaced by h v in the 
entire proof. To do this, we first prove the strong two-scale convergence of u n (x) to u°(x, y) in 



fio, the strong convergence in f2\fi being standard (see footnote [5] on page 30). By Lemma 
the extensions vP e of the restrictions of u v to QJl are bounded in H l (Qo), and therefore we 
can restrict T' to a subsequence that converges strongly to a function w(x) in L 2 (Qo). By 
Theorem 1.14 (a) of [1], we can restrict T' further so that u2(x) also two-scale converges to 
a function in L 2 (f2 ), namely u®(x). These two facts imply the weak L 2 convergence of ul(x) 
to both w(x) and u®{x), from which we obtain 



lim 



In a similar manner, we obtain 



{vP e {x)) 2 dA{x) 



Q.0 



<(x) 



un 



(167) 



lim 



uV{x)YdA{x) 



n 
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Now that we have the strong two-scale convergence of u2(x) to u®(x) and ul(x) to u®(x) 
as well as the two-scale convergence of u^(x)x e (x /rj) to u®(x)x c (y) and ttf (x)xi(x/r]) to 
u®(x)xi(y), we apply Theorem 1.8 of [1] (see footnote^ on page 30) to obtain 



lim 



n 



77^0 



u v {x)) dA(x) = lim / ({^(a;)) Xc( x /v) dA(x) + lim / Xi( x /v) dA(x 



n 



>e(^)) 2 / Xc(y)dA( 2/ )dA(x)+ / K(x)) 2 / Xi (y)dA(y)(L4(x) 



n 



(u° e (x,y)) 2 dA(y)dA(x). (169) 



By definition of m* 7 , if m° < 1 then H^H^o) = 1 for 77 sufficiently small, and by the strong 
two-scale convergence of vP to u°, we obtain 



\u 



1 (m < 1). 



(170) 



Thus u c 7^ because of the relation (164). In addition, u e solves the homogenized system 



(67) with forcing m°f. If in fact, if this solution is unique, we conclude that m° 7^ 0. We will 



show momentarily that it is unique at least in the case that lm(fx*(x)) > 0. We therefore 
obtain 

hP = (m^V (m )~V =: h°, (171) 

obtaining a uniform L 2 bound on h v . All of the results in this proof are therefore valid if 
we replace all occurrences of u with h. Because of the uniqueness of the solution to the 
two-scale variational problem, we conclude the convergences obtained in this proof hold for 
the entire sequence T. 



Finally, we explain why the solution to the homogenized system (153) is unique if 



lm(fi*(x)) > 0. The argument is standard, so we will be brief. Setting the right-hand 



side of the second equation in (153) to zero and using v = h® as a test function, we obtain 

L V ± /i°(x)-V ± he(x) - u 2 fi*(x)\h c (x)\ 2 ] dA(x) = 0. (172) 



e ix 



By its definition (75), Ime*(x) > also, and it follows that h®(x) = 0. 



Theorem 8 The sequence h v (x) converges strongly in the spaces V v to its two-scale limit 
h°(x,y) + r]h}(x,y) in the sense that 

Km\\W(x)-h (x,x/ V )\\ L2{n) = (173) 

and 

lim \\Vhi(x) - V x h°{x,x/rj) - W y h\x,x/r])\\ L2(Q) = 0. (174) 
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Proof. The first limit follows from the strong two-scale convergence of h v (x) to h°(x,y) 



(169) and Theorem 1.8 of pQ. To prove the second limit, we write the expansion 



b (h v (x) — h°(x, y) — r/u 1 ^, y), h v (x) — h°(x, y) — rjh 1 ^, y)) 

= b {h?(x), hP{x)) - b (hP{x), h°(x, y) - r)h x {x, y)) + 
- b (h°(x, y) - rfh 1 ^, y), h v (x)) + b (h°(x, y) - r)u\x, y), h°(x, y) - rfh l (x, y)) . (175) 

The first term on the right-hand side is equal to u 2 c (h v (x), h v (x)) + / (h v (x)), which, by the 
strong two-scale convergence of h v {x) to h°(x,y) tends to 



u 2 / if(x) (h° c (x)) 2 dA(x) + f {h° e (x)) 



(176) 



By using the admissibility of h° + rjh 1 as a test function for two-scale convergence together 
with Theorem 1.8 of pQ and its analog for two-scale convergence on hypersurfaces, we can 
pass to the two-scale limit in each of the other three terms, which is (plus or minus depending 
on the sign in ( 175[ )) 



/ / e(x,y)- 1 [V^h°(x,y)+V^h 1 (x,y)]-\v^h%x,y)+Vj-V(x,y)\ dA(y)dA(x) + 
Jo, Jo L J 



— IU 



a(x,y) 1 ds(y) 



r)G 



\h° e (x) -h?(x)\ 2 dA(x). (177) 



By Theorem 7 and the definition of /i*, this is equal to (176), and we find that the right-hand 



side of ( 175| tends to zero with rj. By the coercivity of b (98), the limit (174) now follows. 



Theorem 9 The reflection and transmission coefficients for the problems of scattering by 
the micro -structured slabs (Problem^ converge to those of the problem of scattering by the 
homogenized slab (equation 153). 



Proof. This follows from the convergence of the Fourier coefficients of the propagating 
harmonics of the solutions h v to those of the solution h® 



lim / h v (x)e-' l( - m+K)x2 dx 2 = [ h° c (x)e- i{m+K)x2 dx 



(178) 



which is valid because of the strong convergence of the solutions and their gradients outside 
the slab. ■ 
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